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ABSTRACT 



Ultramagnetized neutron stars or magnetars have been invoked to explain several 
astrophysical phenomena. We examine how the magnetic field of a magnetar will decay 
over time and how this decay affects the cooling of the object. We find that for sufficiently 
Q\ • strong nascent fields, field decay alters the cooling evolution significantly relative to similarly 

magnetized neutron stars with constant fields. As a result, old magnetars can be expected to be 
bright in the soft X-ray band. The soft X-ray source RXJ 0720.4—3125 may well be the nearest 
^ • such old magnetar. 
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00 . 1. Introduction 

r"| . The birth parameters - the rotation period (P), the strength of the magnetic field (B) and the 

Onj velocity - of neutron stars are a matter of fundamental interest in our study of neutron stars. Another 

equally fundamental issue is the subsequent evolution of the magnetic field. About a decade ago, it was 
thought that most neutron stars are born with parameters similar to that of Crab, initial period of tens of 



milliseconds (perhaps to hundreds) and magnetic field strength confined to a narrow range of 10 G 



(e.g. Narayan 1987) 



Perhaps the single most important parameter is the strength of the magnetic field, B. This strength 



determines the luminosity, the lifetime and even the nature of the energy loss from isolated neutron stars. 



All these three issues affect observability and hence can be ignored at the risk of severe selection. For 



example, Usov & Melrose (1996) argue that when B exceeds 0.1B cr the pulsar wind may be dominated by 
bound pairs rather than freely streaming electrons and positrons; here B cr = m 2 c? jeh w 4.4 x 10 13 G. Such 
strongly magnetized neutron stars may thus be invisible in the radio sky. 

Over the past decade there has been increasing recognition that there is considerable diversity in the 
magnetic field strengths of neutron stars. In particular, there is circumstantial evidence suggesting the 



existence of "magnetars" - neutron stars with magnetic field strengths in excess of B cr (Thompson fc 



Duncan 19951 ). The evidence primarily comes from observations of soft gamma-ray repeaters (SGRs) and 



long-period pulsars in supernova remnants, especially young supernova remnants ( Vasisht & Gotthelf 1997) 



Kulkarni & Frail (1993) have argued that SGRs are young neutron stars, albeit with some special 



parameters, based on their associations with supernova remnants. They and Kouveliotou et al. (1994) 
estimate the birth rate of SGRs is 10% that of ordinary pulsars. 



In this Letter we discuss the possibility of detection of nearby old magnetars. Such soft X-ray emission 
is expected on various theoretical grounds: Thompson & Duncan (1996) through magnetic-field decay and 
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Heyl & Hcrnquist (1997b) by neutron-star cooling. Here, we will examine the consequences of magnetic-field 



decay in magnetars in the context of detailed models of magnetized neutron-star envelopes (Heyl & 



Hernquist 1998a). The essential idea is as follows: the primary store of energy in a magnetar is the energy 
stored in the magnetic field, Eg ~ B 2 R 3 ; here R is the radius of the neutron star. Assuming that the field 
decays on a time scale, td we derive a luminosity Eb/td- Some portion of this could heat up the surface 
leading to potentially detectable soft X-ray emission. 

The organization of the paper is as follows. In section || we discuss a candidate magnetar, 
RX J0720.4 — 3125 ( Haberl et al. 1997 ). The importance of this object is that it is nearby, a few hundred 
parsecs, and thus must be derived from a substantial population (Kulkarni & van Kerkwijk 1998; hereafter 
KV). This is a crucial point in that it provides primary observational motivation to this paper. In 
section |3| we discuss potential field decay mechanisms and i n section ^ we calculate the expected thermal 
emission in old magnetars. In contrast to earlier work {e.g. Thompson & Duncan 1996), these calculations 
take account explicitly the modification of thermal conductivity by strong magnetic fields (B > 10 11 G) 
( Hcrnquist 1984 , Heyl fc Hernquist 1998a ) and its effect on the structure of the neutron-star envelope ( Heyl 
fe Hcrnquist 1997b, Heyl & Hcrnquist 1998a). Finally in section || we present estimates for the number 
of nearby magnetars that are detectable solely based on their thermal emission powered by magnetic field 
decay. 



2. Soft X-ray emission from Magnetar Candidates 



□ 



Thompson & Duncan (1996) argue that the soft X-ray counterpart of SGR 0526—66 (Rothschild 



Kulkarni & Lingcnfelter 1994) is powered by magnetic field decay. It is in this connection that the nearby 
soft X-ray source, RX J0720.4 — 3125 has drawn considerable attention. The source - a bright, soft X-ray 
(kT = 80 eV) pulsar, P = 8.39 s - was identified in the ROSAT All Sky-Survey by |Haberl et al. (1997) 



Motch fc Haberl (1998) and KV searched for an optical counterpart for this source. KV identify a plausible 



faint blue counterpart. In any case, Motch & Haberl and KV argue that the optical counterpart is so faint 
that RX J0720.4 — 3125 is most likely a nearby (100 pc) isolated neutron star. 



KV and Heyl fc Hernquist (1998b) have discussed in detail the possibility that RX J0720.4 
nearby magnetar. Its soft X-ray emission may be powered by the decay of its magnetic field. 



3125 is a 



3. Avenues for Magnetic-Field Decay 

i 1 

Goldreich & Reisenegger (1992) studied several avenues for magnetic field decay in isolated neutron 

stars: Ohmic decay, ambipolar diffusion and Hall drift. Depending on the strength of the magnetic field, 
each of these processes may dominate the evolution: 
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teall ~ 5 x 10 8 ^ yr (4) 

-t>12 \ Pane / 

where L 5 is a characteristic length scale of the flux loops through the outer core in units of 10 5 cm, T 8 is the 
core temperature in units of 10 8 K and B12 is the field strength in units of 10 12 G. Ohmic decay dominates 
in weakly magnetized neutron stars (B <; 10 11 G), fields of intermediate strength decay {B ~ 10 12 — 10 13 G) 
via Hall drift, and intense fields (B ;> 10 14 G) are mostly strongly affected by ambipolar diffusion. 



Thompson fc Duncan (1996) examine the possibility of obtaining an equilibrium between neutrino cooling 



and heating through magnetic field decay for B ~ 10 15 G and T ~ 10 s K. 

Ambipolar diffusion may be decomposed into a solenoidal component which does not upset the 
equilibrium among the neutrons, protons and electrons and an irrotational component. The field evolution 
through the irrotational component is accompanied by modified (or direct if allowed) URCA reactions 
to restore the equilibrium; therefore, it generally operates more slowly than the solenoidal component. 
It is difficult to construct topologically a solution for the ambipolar diffusion of a decaying field without 
both components; therefore, the rate of ambipolar diffusion is set by the slower of the two processes. A 
strong magnetic field may speed the irrotational component by catalyzing the direct URCA process in the 



out-of-equilibrium component (Lcinson & Perez 1997); therefore, we treat both the case where the fast 



solenoidal mode and where the slower irrotational mode set the field decay timescale. 



Calculations 



To estimate the total energy released in the crust by the decay of the magnetic field, we assume that 
the field is a dipole outside of the neutron star and in the insulating envelope, and that the neutron star 
crust is thin relative to the radius of the star. The total energy contained in the neutron star's magnetic 
field is estimated by 

E B = ^B 2 p R 3 (5) 



where B p is the strength of the magnetic field at the pole at the surface, and R is the radius of the neutron 
star. Substantial magnetic energy may be contained in small structures below the envelope; consequently, 
this is a lower bound on the total energy of the neutron star's magnetic field. The power released by the 
decay of the field is 



E r 



(6) 



We will assume that the magnetic energy is dissipated below that region of the crust which supplies the 
bulk of the opacity (i.e. the sensitivity region); the field decay heats the isothermal core of the neutron star. 

To calculate the thermal evolution of the neutron star, we use the T c — T g relationship presented in 



Heyl & Hernquist (1997b) using the conductivities derived by Hernquist (1984), and extended to weaker 



fields using the numerical techniques outlined by Heyl & Hernquist (1998a). The envelope calculations do 



not include a treatment of the Coulomb pressure terms in the equation of state (e.g. Van Riper 1988). 
Although this may introduce uncertainties to the envelope calculations for low effective temperatures 
(T c ff < 10 6 K), the correction to the T c — T e g relationship due to Coulomb interactions is probably not 



a large as implied by the results of Van Riper (1988). The modified Coulomb prescription used by Van 



Riper (1988) yield a large correction to the T c — T e ff relation for precisely those models which contain a 



region of negative total pressure. 



-4- 



We use the same cooling model as Heyl & Hernquist (1997b) coupled with an equation for the field 
decay of the form 

dt P (tohmic ^ambip ^Hall ^ ^ ^ 

The middle term dominates the sum except during the first ~ 500 yr of the decay of the 10 14 G field 
when Hall drift has a comparably strong effect. We study three initial field strengths B Pt o — 10 14 , 10 15 
and 10 16 G, and begin the integrations at t — 1 yr and a core temperature of 10 9 K until t = 10 7 yr using 
logarithmically spaced timesteps. 

Figure [l] shows the evolution of the effective temperature and magnetic field strength for the various 
initial field strengths and assumptions regarding the disposition of the magnetic energy. The release of 
magnetic energy heats the core, increasing the emergent flux from the neutron star. 

For the solenoidal mode, the effect is substantial for an initial field of 10 16 G at t = 10 3 yr, the effective 
temperature is 40% higher than if the field did not decay. For Bq — 10 15 G, the effective temperature is 10% 
higher at t = 10 4 yr. Because the initial magnetic field strength has not been assumed to be infinite, the 
evolution of the magnetic field strength departs from a power law. However, for the strongest initial field 
(Bq = 10 16 G), the evolution approaches a power law after a few characteristic decay timescales. The slope 
of the power law changes appreciably when the neutron star leaves the neutrino-dominated cooling era, and 
photon emission begins to dominate the cooling (at t ~ 10 6 yr). Additionally, as the energy released in the 
field decay heats the core, the decay is slowed slightly as apparent in the upper-right panel of Figure [l| 

The character of field decay in the irrotational mode is different. Since the field tends to decay 
more slowly, its effect on the cooling of the neutron star is most pronounced later in the star's evolution 
(t ~ 10 6 yr). For Bq = 10 16 G, the neutron star with a decaying field is brighter throughout its evolution 
than a similar neutron star with a static field. For a more weakly magnetized neutron star (Bq = 10 15 G), 
we find that heating by the decaying field and neutrino cooling attain a near equilibrium for t ~ 10 5 yr, 
similar to the result found by Thompson fc Duncan (1996)| . Regardless of their initial fields as long as 



Bq ;> 10 15 G, neutron stars with decaying fields are much brighter during the photon-cooling epoch than 
their coevals with static fields. 

For weaker fields however (B <J 10 14 G), both the total energy contained in the magnetic field and 
the rate of its decay are insufficient to perturb the cooling evolution during this epoch. Even if ambipolar 
diffusion is unimportant, the field decay by Hall drift can result in a significant increase in the surface 
emission. 

Our models for the decay and radiation of the magnetic field energy yield comparable X-ray luminosities 
to the calculations of Thompson fc Duncan (1996)| who discuss the mechanisms for field decay in magnetars 



in further detail but extrapolate the envelope models Van Riper (1988 j| which were derived for more weakly 



magnetized objects (B < 10 14 G) to B ~ 4.4 x 10 15 G. For field decay dominated by the solenoidal mode, we 
generally find that photon cooling begins dominate earlier than for irrotational decay. For the irrotational 



mode our results compare well with Thompson fc Duncan (1996) conclusions 



5. Discussion 



In this paper we have estimated the cooling radiation from magnetars. Two factors conspire to make 
them brighter relative to the standard pulsars. First, the decay of the magnetic field by ambipolar diffusion 
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Fig. 1. — The upper panels depict the evolution of the effective temperature (left) and magnetic field (right) 
for the case where the magnetic field decays through the solenoidal mode. The lower panels depict the 
results for decay through the irrotational mode. The bold solid, dotted and dashed lines give the results for 
Bi = 10 14 , 10 15 and 10 16 G respectively. The light lines give the results without any field decay, and the 
long-short-dashed line gives the cooling evolution for an unmagnetized neutron star. 
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(see section 0) is a new source of heat. Second, thermal conduction is affected by strong magnetic fields 



(Heyl & Hernquist 1997b). The main uncertainties are the rate at which the field can decay which depends 
on how the neutron-electron-proton plasma returns to equilibrium. If the solenoidal mode operates, neutron 
stars with decaying fields of 10 15 to 10 16 G may be 50% to 250% more luminous than objects with static 
fields similar in magnitude during the neutrino-cooling epoch. If the irrotational mode dominates the 
cooling, the effect is strongest later in the neutron stars evolution, during the photon-cooling epoch. 

Here, we have considered only envelopes consisting of iron. If neutron star envelope consists largely 
of lighter elements (hydrogen or helium), the opacity of the envelope is even smaller, making the putative 



magnetars even brighter (Heyl & Hernquist 1997a, Potekhin, Chabrier & Yakovlev 1997). By increasing 
the transparency of the envelope (e.g. Heyl & Hernquist 1997a), even a static field of B ~ 10 15 (weaker 



than required by Thompson & Duncan (1996)) can account for the quiescent emission of SGRs. Regardless 
of composition though, neutron stars with decaying magnetic fields will be more luminous than their 
non-decaying coevals. 

We take RX J0720.4— 3125 as a more modest exemplar (as compared to the SGRs) for magnetar class; 
it is a remarkably bright source in the ROSAT All-Sky Survey with a PSPC count rate of 1.67 s _1 . The 
distance to this object is uncertain but it is estimated to lie between 100 to 400 pc. If we take a nominal 
limiting flux for the ROSAT All-Sky Survey Bright Source Catalogue of 0.10 s~ x where it has a sky coverage 
of 92%, we would expect ~ 50 similar objects to be appear in the survey. However, the observed population 
will be restricted severely by the opacity (photoelectric absorption) of the neutral interstellar medium. 
Nonetheless, it is not unreasonable to expect a dozen or so of objects similar to RX J0720.4 — 3125 but 
fainter. 

For the reasons discussed above, in a flux-limited and a distance-limited sample of soft X-ray sources, 
middle-aged (~ 10 yr) magnetars will be over-represented relative to middle-aged neutron stars with typical 
field strengths of 10 12 G (or weaker). In principle, it is possible to identify this population. However, in 
practice, it is difficult to identify weaker cousins of RX J0720.4 — 3125 using the existing data. Fortunately, 
AXAF observations of the ROSAT sample of soft X-ray sources will have the necessary astrometric precision 
and spectroscopic resolution differentiate cooling neutron stars from dominant contaminating populations 
of white dwarfs, soft AGN and CVs. 

Magnetic field decay along with enhanced thermal conduction makes ultramagnetized neutron 
stars much brighter than their weakly magnetized coevals during the neutrino-dominated cooling era; 
consequently, the search volume for young neutron stars associated with supernova remnants is much larger 
for magnetars with decaying fields than for weakly magnetized neutron stars. These unusual objects may be 
commonplace in our flux-limited view of the galaxy, and these processes may explain the quiescent emission 
of anomalous X-ray pulsars and soft-gamma repeaters among other phenomena. 



The authors would like to thank P. Goldreich for useful discussions. J.S.H. acknowledges a Lee A. 
DuBridge postdoctoral scholarship. 
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